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ICE  ARCHING  AND  THE  DRIFT  OF  PACK  ICE 
THROUGH  RESTRICTED  CHANNELS 


D.S.  Sodhi 


INTRODUCTION 

In  this  report  the  factors  influencing  ice  arching 
across  straits  or  gulfs  are  investigated  by  developing 
an  analogy  between  the  drift  of  pack  ice  through  a 
restricted  channel  between  converging  land  masses 
and  the  gravity  flow  of  granular  material  through 
hoppers  or  chutes.  This  approach  appears  reasonable 
because  of  the  similarity  of  the  basic  mechanics  of 
the  two  situations,  pack  ice  being  composed  of  inter- 
acting ice  floes  that  transfer  stresses  across  their 
boundaries  in  a manner  that  is  presumably  similar  to 
the  force  interaction  between  grains  of  a granular 
material.  Also  supporting  the  applicability  of  such  an 
approach  is  the  knowledge  that  Janssen’s  (1895)  theory 
for  the  flow  of  grain  in  a silo  has  been  used  (Pariset  and 
Hauser  1961,  Michel  1968)  to  calculate,  for  fragmented 
ice  covers  moving  through  straight  river  channels, 
forces  that  were  in  reasonable  agreement  with  experi- 
mental observations. 

Other  work  of  interest  is  by  Shapiro  and  Burns 
(1975)  who  described,  via  the  use  of  LANDS  AT 
imagery,  a short-term  episode  (6-8  March  1973)  of 
the  rapid  deformation  of  pack  ice  in  the  southeastern 
Chukchi  Sea  and  its  subsequent  transport  south  through 
the  Bering  Strait.  They  suggested  that  the  process  of 
ice  passing  through  the  strait  is  mechanically  similar 
to  the  extrusion  of  metal  through  a restricted  orifice. 
They  did  not,  however,  attempt  to  use  extrusion  flow 
theory  to  make  detailed  calculations  concerning  the 
behavior  of  pack  ice.  We  believe  that  granular  media 
models  have  advantages  here,  inasmuch  as  pack  ice  is  a 
granular  medium  both  before  and  after  passing  through 
a strait  while  a metal  is  invariably  continuous  after 
passing  through  a die.  Also,  the  principal  force  which 
drives  the  ice,  the  force  of  the  wind,  can  be  considered 
more  analogous  to  the  specific  weight  of  a granular 
material  than  to  the  forces  involved  in  extrusion. 


A kinematic  model  study  of  the  arching  phenomena 
in  drifting  ice  was  undertakn  by  Calkins  and  Ashton 
(1975)  who  performed  experiments  in  a flume  using 
square  plastic  pieces  to  simulate  ice  floes.  The  plastic 
floes  were  released  at  a controlled  rate  upstream  from 
a channel  restriction  and  the  occurrence  or  nonoccur- 
rence of  a stable  arch  was  noted.  Good  correlations 
were  obtained  between  the  threshold  of  arching  and 
parameters  such  as  the  supply  rate  of  the  surface  area 
of  ice,  the  surface  discharge  at  the  gap,  and  the  ratio 
between  the  size  of  the  individual  floes  and  the  size  of 
the  gap.  Because  their  interests  were  directed  toward 
ice  jams  in  rivers,  Calkins  and  Ashton  (1975)  did  not 
study  the  case  when  the  ratio  of  the  floe  size  to  the 
gap  dimension  is  small,  the  subject  of  interest  here. 

In  this  report  a theory  will  be  developed  and  then 
applied  to  ice  arching  in  the  Bering  Strait  and  Amundsen 
Gulf  via  an  analysis  of  the  geometry  of  the  arching  phe- 
nomena as  revealed  by  satellite  imagery.  This  analysis 
allows  one  to  calculate  representative  values  for  the 
angle  of  internal  friction  and  the  cohesive  strength  of 
pack  ice  as  viewed  on  a large  scale.  In  addition,  com- 
parisons are  made  between  the  lead  patterns  as  revealed 
by  the  satellite  imagery  and  faults  (tangential  velocity 
discontinuities)  that  are  determined  by  a radial  stress 
field  solution  of  the  basic  equations. 

THEORY 

An  essential  requirement  of  the  various  theories  for 
the  flow  of  granular  material  is  that  they  obey  Coulomb's 
yield  criterion.  One  group  of  theories  (Drucker  and 
Prager  1952)  is  based  on  the  assumptions  that  1 ) the 
Mohr-Coulomb  yield  function  is  the  plastic  potential 
for  the  material,  so  that  the  strain  rate  components  are 
related  to  the  yield  function  by  the  associated  flow  rate 
of  plasticity,  and  that  2)  the  principal  axes  of  the  stress 
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and  strain  rate  tensors  coincide.  An  inevitable  conse- 
quence of  these  postulates  is  that  there  is  an  unrealistic 
volume  expansion  during  deformation.  In  an  attempt 
to  overcome  this  difficulty,  work  hardening  models 
have  been  proposed  by  Drucker  et  al.  (1957), 

Jenike  and  Shield  (1959),  and  jenike  (1961 ).  In 
these  models,  cohesion  is  a function  of  the  density  of 
the  material  as  well  as  the  hydrostatic  pressure  to 
which  it  is  subjected.  Another  group  of  theories 
(de  Josselin  de  Jong  1959,  Spencer  1964,  Mandl  and 
Lugue  1970)  has  sought  to  retain  volume  constancy, 
thus  abandoning  the  concept  of  normality  and  the 
coincidence  of  the  principal  axes  of  the  stress  and 
strain  rate  tensors.  Of  these  latter  theories,  that  of 
Spencer  (1964)  appears  to  be  the  most  satisfactory, 
particularly  since  it  allows  faults  on  the  Coulomb 
surfaces,  and  has  recently  been  extended  to  include 
gravitational  and  inertial  acceleration  terms  (Morrison 
and  Richmond  1976). 

No  matter  which  theory  is  followed,  one  ends  up 
with  a set  of  hyperbolic  differential  equations  which 
can  be  solved  by  the  method  of  characteristics  (slip 
line  fields),  provided  the  boundary  conditions  are 
well-defined.  Simplifying  assumptions  commonly 
have  to  be  made  to  arrive  at  simple  solutions  and  a 
number  of  such  solutions  exist  in  the  literature 
(Walker  1966,  Johanson  and  Jenike  1967,  Lee  1973). 

In  this  report,  we  would  like  to  pay  attention  to 
three  aspects  of  the  problem:  1 ) the  limiting  span  of 
an  arch  that  has  formed  across  a strait,  2)  the  flow  of 
pack  ice  through  converging  channels,  and  3)  the 
critical  condition  of  stoppage  of  such  a flow.  The 
theories  dealing  with  each  aspect  of  the  problem  are 
discussed  briefly  below. 

Limiting  span  of  an  arch 

The  theory  for  predicting  the  limiting  spans  for  the 
arching  of  granular  material  is  based  on  a static  analysis 
of  the  stress  distribution  within  the  material  to  ascertain 
whether  the  state  of  stress  is  always  below  the  critical 
state.  To  derive  expressions  determining  the  upper 
and  lower  boundaries  for  the  limiting  span  over  which 
a free  arch  can  form,  we  will  follow  the  analysis  given 
by  Richmond  and  Gardner  (1962)  and  Gardner  (1962) 
by  making  the  following  modifying  assumptions: 

1 . In  the  vicinity  of  a free  arch  the  stress  is  invariate 
with  respect  to  the  longitudinal  coordinate.  This  as- 
sumption has  been  verified  by  many  investigators  in 
deep  vertical  grain  silos  and  by  Michel  (1968, 1971)  in 
long  straight  channels  full  of  fragmented  ice. 

2.  The  material  can  be  approximated  as  a continuum 
which  obeys  Coulomb’s  failure  criterion. 


3.  The  failure  criterion  holding  at  the  boundary 
between  the  granular  material  and  its  container  (i.e. 
between  the  pack  ice  and  the  “containing”  fast  ice) 
is  the  same  as  within  the  material. 

In  the  following,  the  upper  boundary  will  be  defined 
as  the  maximum  span  across  which  an  arch  can  possibly 
form  with  the  state  of  stress  at  all  points  still  remaining 
below  the  critical  level.  The  lower  boundary,  on  the 
other  hand,  is  the  minimum  span  through  which  flow 
can  still  occur. 

To  develop  the  conditions  requisite  for  free  arching, 
we  will  choose  the  Cartesian  coordinate  system  shown 
in  Figure  1 . Making  use  of  assumption  1 , the  equilibri- 
um equations  can  be  satisfied  by  the  following  ex- 
pressions for  the  normal  and  shear  stress  components 
ox  and  rxy: 

ox=-p  (1) 

and 

tx  y = T*  (2) 

where  p is  a constant  and  y is  the  body  force  per  unit 
area.  The  state  of  stress  at  a point  on  the  free  arch  can 
be  conveniently  expressed  with  the  help  of  a Mohr 
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circle  (Fig.  2)  which  must  pass  through  the  origin  of 
the  (r,  o)  graph.  Coulomb’s  failure  criterion  is  ex- 
pressed by  a yield  line: 

t = c + o tan  <(>  (3) 

where  c and  ft  are  the  cohesive  strength  and  angle  of 
internal  friction,  respectively.  If  the  Mohr  stress  circle 
is  tangent  to  the  yield  line,  the  state  of  stress  is  critical. 
According  to  eq  2,  the  maximum  shear  stress  occurs  at 
the  boundary,  and  we  can  express  the  upper  and  lower 
bounds  by  points  U and  L in  Figure  2 by  incorporating 
assumptions  2 and  3.  The  limiting  values  of  the  span 
are  given  by  the  following  expressions  (Richmond  and 
Gardner  1962): 


Figure  3.  Effective  yield  locus  in 
jenike 's  (1964)  model. 
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The  shape  of  the  limiting  arch  can  be  obtained  with 
the  aid  of  the  stress  circle  shown  in  Figure  2.  If  6 is 
the  angle  between  the  horizontal  and  the  tangent  to 
the  free  surface,  the  arch  profile  is  given  by 


= 32L.  (6) 

dx  ox  p 

Integrating  and  assumingy  = 0 when  x = we  obtain 

'■#['-(*)']  171 

indicating  that  the  profile  of  the  limiting  arch  should 
be  a parabola. 

The  angle  0b  in  Figure  1 will  either  be  equal  to  0g 
or  0L  in  Figure  2,  depending  upon  whether  the  arch 
is  of  the  upper  or  lower  boundary  type.  0y  is  equal 
to  45°  and  0L  can  be  shown,  with  the  help  of  the 
stress  circle  in  Figure  2,  to  be  equal  to  rr/4  +■  0/2.  The 
base  angle  0b  can  therefore  be  used  to  estimate  the 
angle  of  internal  friction  (ft. 

Flow  of  pack  ice  through  converging  channels 
When  ice  is  moving  through  a restricted  channel, 
Jenike’s  (1964)  model  of  the  steady  gravity  flow  of 
granular  material  through  a converging  channel  can  be 
used  to  obtain  the  locations  at  which  velocity  dis- 
continuities (i.e.  cracks  or  leads)  will  occur.  In  this 


model  there  is  assumed  to  exist  a family  of  similar  but 
varying  yield  loci  corresponding  to  different  degrees 
of  consolidation.  Each  locus  terminates  at  the  stress 
appropriate  to  its  consolidation  (i.e.  the  material  can- 
not be  subjected  to  a higher  stress  without  becoming 
further  consolidated).  The  envelope  of  terminal  Mohr 
circles  of  individual  yield  loci  that  is  termed  the 
“effective  yield  locus”  (Fig.  3)  represents  the  shear- 
compressive  stress  characteristic  of  the  material.  For 
many  granular  materials,  the  effective  yield  locus  is 
linear  and  passes  through  the  origin  (Jenike  1961 ); 
therefore,  we  will  also  assume  that  this  is  true  for  the 
case  of  pack  ice.  When  the  state  of  stress  is  such  that 
the  Mohr  circle  is  tangent  to  the  effective  yield  locus, 
the  material  is  said  to  be  in  a critical  state,  i.e.  steady 
state  of  flow.  According  to  this  theory,  one  obtains 
two  separate  sets  of  stress  and  velocity  characteristics, 
and  any  velocity  discontinuity  will  be  along  the  velocity 
characteristics.  In  the  application  considered  later,  we 
will  use  the  radial  stress  field  solution  proposed  by 
Jenike  (1964),  which  is  valid  near  the  outlet  of  the 
converging  channel,  and  assume  that  the  hydrostatic 
stress  is  proportional  to  the  distance  from  the  vertex 
of  the  channel. 

Stoppage  of  flow 

The  most  complete  analysis  of  the  non-steady  flow 
of  granular  material  in  a converging  channel  is  by 
Morrison  and  Richmond  (1976).  It  uses  the  ideal  soil 
model  proposed  by  Spencer  (1964),  which  is  based  on 
the  assumptions  of  material  isotropy  and  incompres- 
sibility together  with  the  postulate  that  deformation 
occurs  by  shear  along  the  characteristic  curves  of  the 
stress  equations.  This  model  is  particularly  satisfactory 
for  the  present  application  because  it  allows  tangential 
velocity  discontinuities  (faults)  along  the  stress 
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Figure  4.  Converging  channel. 


characteristics.  Based  on  Morrison  and  Richmond’s 
(1976)  analysis,  the  critical  value  of  the  slope  B of  a 
converging  channel  (see  Fig.  4)  for  which  no  flow  or 
arching  may  occur  is  given  by 


fi  -f  cqt0  (1  - 3 sin0) 
c (pg)  WQ  (l  - sin  0) 


where  fic  = critical  value  of  the  slope  of  the  channel 
c = cohesive  strength  of  the  material 
0 = angle  of  internal  friction  of  the  material 
W0,  W,  = half  widths  of  channel  at  entry  and  at  exit, 
respectively 
R = W0/W, 

(pg)  = specific  weight  of  tTie  material 
m = (2  sin  0)/(1  - sin  0). 

If  B < Bc  arching  will  occur. 


APPLICATIONS 
St.  Lawrence  Island 

This  island  is  situated  directly  south  of  the  Bering 
Strait.  As  the  ice  drifts  southward  through  the  strait 
and  moves  past  the  island,  there  is  a wedge-shaped 
accumulation  of  ice  that  builds  up  on  the  north  side  of 


St.  Lawrence  Island 


Figure  5.  LANDS  A T imagery  showing 
wedge-shaped  ice  accumulation  north 
of  St.  Lawrence  Island,  24  February 
1975. 


the  island,  as  shown  in  Figure  5.  This  phenomenon  is 
similar  to  the  piling  of  granular  material  on  an  ob- 
struction placed  in  the  flow  of  granular  material. 

Using  Spencer’s  (1964)  soil  theory,  the  sliding  sur- 
faces can  be  interpreted  as  the  stress  characteristic  lines 
for  a constant  orientation  of  principal  stresses  in  the 
pack  ice,  as  shown  in  Figure  5.  The  smaller  angle  be- 
tween the  characteristic  lines  can  be  shown  to  be 
equal  to  [(rr/2)  - 0] . The  wedge  angle  in  Figure  5 is 
approximately  equal  to  59°,  trom  which  we  can  cal- 
culate the  angle  of  internal  friction  to  be  equal  to  31°. 


A similar  wedge-shaped  ice  buildup  has  been  noted 
by  Dunbar  and  Weeks  (1975)  in  the  ice  next  to  Prince 
Edward  Island  in  the  Gulf  of  St.  Lawrence.  The  angle 
between  the  characteristic  lines  at  that  site  was  50°, 
corresponding  to  an  angle  of  internal  friction  of  40°. 


Amundsen  Gulf 

The  ice  arch  across  Amundsen  Gulf  in  the  western 
Canadian  Arctic  (»  71  °N,  123°W)  is  quite  clearly  shown 
in  the  LANDSAT  imagery  obtained  on  12  and  1 3 June 
1973  (Fig.  6 and  7).  The  profile  of  die  arch  is  approxi- 
mately parabolic,  which  is  in  agreement  with  the  shape 
predicted  by  eq  6.  Although  the  breakup  of  the  arch 
had  started  on  12  June  (Fig.  6),  the  complete  collapse 
took  place  on  13  June  (Fig.  7).  The  shearing  of  ice  at 
the  base  of  the  arch  is  evident  in  Figures  6 and  7. 

The  angle  of  internal  friction  0 is  obtained  by 
measuring  the  angle  of  the  arch  at  the  base,  which  is 
equal  to  2/t/X  in  the  case  of  a parabola,  h being  the 
height  and  X being  the  semi-span  of  the  parabola. 
Because  this  angle  is  greater  than  45°,  it  suggests  that 
the  arch  is  a limiting  lower  bound  arch.  This  is,  of 
course,  in  agreement  with  the  observation  that  the 
arch  had  started  to  break  up.  From  the  imagery,  the 
angle  of  the  arch  at  the  base  is  found  to  be  approxi- 
mately 62.5°,  from  which  0 is  estimated  to  be  2(62.5  - 
7t/4)  = 35°.  Now  that  0 is  known,  a value  can  be  ob- 
tained for  the  cohesive  strength  c of  the  pack  ice  via 
eq  5 if  an  adequate  estimate  can  be  given  of  y,  the 
body  force  per  unit  area.  In  the  sea  ice  case,  y cor- 
responds to  the  resultant  force  on  the  ice  primarily 
produced  by  the  combination  of  the  wind  force  (on 
the  upper  ice  surface)  and  the  water  force  (on  the 
lower  ice  surface).  As  a first  approximation  we  will 
consider  the  water  force  to  be  negligible  and  discuss 
only  the  wind  force. 

The  closest  wind  observations  to  the  area  of  interest 
were  from  Sachs  Harbour,  N.W.T.,  on  12  June  1973 
where  the  maximum  and  mean  winds  were  reported 
to  be  6.25  m/s  and  4.5  m/s,  respectively,  from  the 
southeast  to  the  south.  The  direction  of  the  waves  in 
the  LANDSAT  imagery  of  12  June  indicates  that  the 
wind  direction  was  southerly.  The  direction  of  drift  of 
ice  floes,  as  determined  from  their  positions  on  two 
consecutive  days,  suggests  that  the  average  wind  direc- 
tion on  13  June  1973  was  southeasterly,  approximately 
perpendicular  to  the  base  of  the  parabola.  A rough 
estimate  of  the  wind  stress  per  unit  area,  which  we  will 
take  to  be  equal  to  y,  can  be  obtained  from 


y = pDv 2 


where  p is  the  mass  density  of  the  air,  D the  drag  co- 
efficient, and  v the  wind  velocity.  Using  a drag  coef- 
ficient of  0.001  (Hiblcr  and  Mock  1973,  Banke  and 


Smith  1973),  the  wind  stress  is  estimated  to  be  equal 
to  0.051  N/m2.  Now  taking  the  span  2X  to  be  approxi- 
mately 123  km,  we  can  estimate  from  eq  5 that  the 
cohesive  strength  is  1993  N/m.  It  is  interesting  to  note 
that  Pariset  and  Hauser  (1961)  estimated  the  cohesive 
strength  between  fragmented  river  ice  and  the  river 
bank  to  be  equal  to  a value  between  1095  N/m  and 
1320  N/m,  values  which  are  of  the  same  order  of 
magnitude  as  we  obtained. 


Bering  Strait 

Figures  8 and  9 show  infrared  imagery  of  the  Bering 
Strait  region  obtained  by  the  U.S.  Defense  Meteorology 
Satellite  Program.  In  these  images,  lighter  areas  indicate 
high  temperature  zones,  which  in  this  case  are  either 
open  or  newly  refrozen  leads  in  the  pack  ice.  it  is 
obvious  from  these  images  that  the  Bering  Strait  is,  in 
fact,  a converging  channel  similar  to  those  treated  in 
the  chute  and  hopper  theories.  (Another  view  of  this 
region  is  shown  on  the  cover.) 

Using  the  results  of  Jenike’s  radial  stress  field  solu- 
tion (Johanson  and  Jenike  1967),  which  we  discussed 
earlier,  velocity  characteristics  were  constructed  for 
the  converging  channel  shown  in  Figure  10  and  were 
matched  with  the  major  leads  (cracks)  shown  in  the 
LANDSAT  imagery.  A good  fit  was  obtained  when 
the  effective  angles  of  internal  friction  within  both 
the  pack  ice  and  the  stationary  fast  ice  were  chosen  to 
be  30°. 

We  can  use  eq  8 to  predict  the  conditions  required 
for  no  flow  or  arching  in  the  strait,  modeled  as  a con- 
verging channel  of  slope  B,  by  replacing  the  specific 
weight  (pg)  with  the  applicable  wind  stress.  If  the  wind 
stress  is  less  than  that  required  to  produce  the  critical 
condition  of  flow,  arching  may  occur.  To  obtain  an 
estimate  of  critical  wind  velocity,  values  of  B,  R and 
W0  are  measured  from  Figures  8 and  9.  Then  using 
the  values  of  c and  0 obtained  earlier,  the  wind  velocities 
necessary  to  cause  B = Bc  are  calculated  from  eq  9.  The 
relations  indicate  that  as  long  as  the  wind  velocity  ex- 
ceeds a value  of  5.6  m/s  through  the  Bering  Strait,  there 
will  be  a flow  of  ice  through  the  Strait  in  the  general 
manner  shown  in  the  imagery. 

In  some  of  the  DMSP  images  of  the  Bering  Strait 
region  (Fig.  1 1 ),  one  can  observe  long  leads  in  the  south- 
eastern Chukchi  Sea  that  separate  the  moving  pack  ice 
from  the  stationary  ice.  One  can  also  observe  smaller 
leads  that  originate  from  the  long  leads.  These  are 
oriented  in  a direction  perpendicular  to  the  extensional 
strain  in  a manner  similar  to  the  orientation  of  crevasses 
in  a glacier.  If  we  assume  that  the  pack  ice  deformation 
in  Figure  1 1 occurring  along  the  long  leads  is  basically 
in  pure  shear,  making  these  leads  one  of  the  stress  charac- 
teristics, then  it  can  be  shown  that  the  angle  between 
the  stress  characteristic  and  the  line  perpendicular  to 
the  direction  of  extensional  strain  is  (ff/4)  - (0/2). 
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Figure  / /.  DMSP  imagery  showing  crack  pattern  due  to  shear- 
ing, Bering  Strait,  8 December  1975. 


As  measured  on  the  DMSP  imagery  this  angle  is  found 
to  be  approximately  30°,  corresponding  to  an  angle  of 
internal  friction  0 of  also  30°. 


CONCLUSION 


encouraging  in  that  there  is  correspondence  between 
observed  arching  and  lead  patterns  and  calculated  pat- 
terns. Also,  values  determined  via  the  model  tor  the 
angle  of  internal  friction  and  cohesive  strength  appear 
quite  reasonable  when  compared  with  similar  values 
available  in  the  literature. 


In  the  preceding,  models  originally  developed  to 
describe  the  arching  or  the  movement  of  granular 
materials  such  as  grain  or  sand  through  hoppers  or 
chutes  have  been  applied  to  the  arching  and  drift  of 
pack  ice  on  a scale  of  50  to  500  km.  Inasmuch  as 
experimental  verification  of  such  theories  as  applied  to 
sea  ice  is  not  presently  possible,  comparisons  arc  made 
with  ice  deformation  patterns  as  observed  with  satellite 
imagery.  The  results  of  these  comparisons  are  quite 
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